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The ultraviolet extragalactic background light: dust 
extinction and the evolution of the cosmic star formation 
rate from z = to ~ 0.6 
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ABSTRACT 

We show that the accumulated light of galaxies in the ultraviolet can be evaluated from 
their luminosity density as a function of the evolution of the cosmic star formation 
rate and dust extinction properties. Constraints on the evolution rate are expected in 
future. Data available at the moment are consistent with an evolution rate at low z 
steeper than (1 + z) 3 ' 5 . A shallower rate remains possible if the luminosity-weighted 
dust extinction at 2000 A, as suggested by some data, is lower than ~ 1.2. 
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1 INTRODUCTION 

The present-day luminosity density of galaxies and evolution 
effects are folded into the cumulative emission from galax- 
ies. Disentangling these quantities is known to be difficult 
because of the need to account for the spectral energy distri- 
butions of the galaxies at wavelengths shorter than the win- 
dow of observation. The cumulative emission from galaxies 
may be obtained directly from galaxy number counts or iso- 
lated within gl obal measurements of the diffuse back ground 
radiation (e.g. iBernstein. Freedman fc Made ira 2002). 

In the far-UV wavelength range, the situation is some- 
what more simple. The spectral energy distributions are in 
first approximation dominated by star formation activity, 
avoiding the need to keep track of too many categories of 
galaxies and opening the possibility of constraints on the 
history of stellar birth in galaxies. The Lyman break enters 
the window of observation at relatively low redshift, reduc- 
ing the amount of look-back time involved. The accumulated 
light of galaxies is the dominant source of extragalactic back- 
ground light which is itself a si gnificant contributor to the 
total background radiation (e. gjM^rtm^Ihur witz fc Bowverl 
ll99ll:lArmand. Milliard fc Deharvenslll994l) . ~ 

Three factors have recently made a specific examination 
worth to attempt. First, the luminosity density of galaxies 
is now available at far-UV wavelengths llTrever et al.lll99c1: 
|Sjffhvan^^l ]|2000l). Second, the g alaxy counts in the UV 
iGardner. Brown fc Fergu son 2000) are now deep enough to 
provide a reliable evaluation of the background radiation due 
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to galaxies, without the difficulties of subtracting uncertain 
components from the observations of the total diffuse back- 
ground radiation. Third, t he measurements of t he luminosity 
density of galaxies at Ha jGallego et alJll995l) . can be con- 
verted into a luminosity density of ionizing photons, then 
connected by stellar population models with the luminosity 
density above the Lyman break and used as a constraint on 
the average spectral energy distribution of galaxies in the 
far-UV. 

In this paper, our motivation is to evaluate the far-UV 
integrated light of galaxies with a limited set of parame- 
ters and to explore whether useful constraints may be de- 
rived on the history of star formation. Despite an explosion 
in the quantity of faint galaxy data, there are yet signif- 
icant disagreements in t he determinations of the star for- 
mation rate (SFR) (e.g. iHopkins et al-lEoOll iHogd l2002tl 
that trace the cosmic star formation history. Even in rel- 
ative terms and at low z, the rate of evolution is still a 
matter of debate, with a parameterization ranging from 
(1 + z) 1 - 5 JCowie. Songaila fc Bareerl Il999t |Wilson et al " 



20021) to (1 + z) 4 (e.g. iLillv et alJ Il99a iMadau et a" 



1996 : lMadau. Pozzetti fc Dickinsonll99ctl for measurements 
based on U V rest fram e data, not to speak of more extreme 
values (e.g. Hogg 2002) when measurements at other wave- 
lengths are included. 



2 FORMULATION 



The cosmological radiative transfer equation jPeebles)ll993l) 
for sources with proper specific volume emissivity e(u, z) (in 
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ergs cm 3 Hz 1 s 1 ) gives a mean specific intensity at ob- 
served frequency vp and for an observer at redshift (e.g. 
iBechtoid et aTlll987h . 

where v = vq(\ + z) and the opacity of the intergalactic 
medium is neglected for our application in the non-ionizing 
ultraviolet and at relatively low redshift. The relation be- 
tween the proper length increment and the redshift incre- 
ment is given by 

=- = -^(i + (i + zf + n A )-°- 5 (2) 

az Ho 

where the cosmological parameters and c have their usual 
meanings. In our case the volume emissivity is the ultravi- 
olet emission resulting from the sta r formation act i vity i n 
galaxies. Following the notations of IBechtoid et aD il987t) . 
we assume that it can be written as 

e(u, z) = e(i/ )(l + zf^{z)s{v/u ) (3) 

where e(fo) is the current local luminosity density of galax- 
ies at frequency vo, s(v/vq) represents the spectral shape 
normalized to 1 at the frequency and tp{z) accounts for 
any proper evolution in the luminosity density. We take vo as 
corresponding t o 1595 A , the pivot wa velength of the far-UV 
galaxy counts of lGardner et all i2000l> . The evolution ip(z) is 
currently parameterized as (1 + z) 7 . For the purpose of sim- 
plification and according to models of star-forming galaxies 
(continuous star formation), we assume that s(v/vq) can be 
parameterized as (v/vo) a between our window of observa- 
tion and the Lyman break. 

In these conditions the integrated light from galaxies at 
1595 A is given by 

^ = i^ i595 >r vizirs dz w 

In practice the upper bound in equation (4) will be approx- 
imately limited to the redshift at which the Lyman break 
reaches the pivot wavelen gth of the observatio ns of the in- 
tegrated light of galaxies iGardner et al. 2000). 



l2002HBuat et ~aTl l2002F) lead us to stop our extrapolation to 
~ 1000 A because of the many absorption features that de- 
press the flux before the Lyman break is effectively reached. 
In these conditions, the effective slope a is defined by the 
slope between 2000 A and 1000 A. Accounting for the dif- 
ferential dust extinction, the effective slope can be written 
as, 

_ _„ j, 0.4(Aiooo — Agooo) 
° ~ ' (log(2000) - log(1000)) 

with Aiooo and A2000 being the amount of extinction at 1000 
A and 2000 A respectively. It is finally found to depend 
on A2000 and the ratio of the reddening law at 1000 A and 
2000 A , fc(1000)/fc(2000) 

a = -0.1 + 1.328A 20 oo(fc(1000)/fc(2000) - 1) (5) 

Following the prescriptions of ICalzetti et alj (|200(J) and 
their e xtrapolation shortward of 1200 A bv lLeitherer st alJ 
(20112 ). we have fc(1000)/fc(2000) = 13.88/8.87 = 1.565 

The value of e(1595) is derived fr om the luminosity den - 
sity of galaxies observed at 2000 A bv lSullivan et alJ {2000), 
taken as 0.84 x 10 38 ergs s" 1 A" 1 Mpc" 3 for H = 70 km 
s -1 Mpc -1 . It is then scaled from z — 0.15 to z = using 
the factor 1/(1. 15) 7 , and from 2000 A to 1595 A using the 
slope a determined as a function of ^2000- 

The integrated light from galaxies, calculated from 
equation (4), is displayed in Fig. 1 as a function of ^2000 
for three values of the evolution rate. It has been converted 
into units of photons cm -2 s -1 A -1 ster -1 for comparison 
with the range of values 144 - 19 5 phot ons cm -2 s -1 A -1 
ster -1 reported bv lGardner et al] i2000t) . An upper bound 
of 2 = 0.6 has been used in equation (4) but the integrated 
light of galaxies would be increased by any flux shortward 
of the Lyman break. The increase has been found to be 
less than ~ 10% (even in the case of the evolution rate of 
(1 + z) 4 ) for an upper limit to the Lyman continuum escape 
fraction of ~ 6% in the nearby star-forming galaxies (e.g. 
[H urwitz et al.lll997t iHeckman et alJboOlt iDeharveng et all 
1200 ll) . 



3 EVALUATION OF THE UV BACKGROUND 
RADIATION DUE TO GALAXIES 

We first discuss the two quantities a and e(1595) that, in 
addition to the evolution factor 7, enter the equation (4) for 
the calculation of the background radiation due to galaxies. 

The normalized spectral shape s(v/vo) is written as 
(y/vo) a = (1 + z) a under the assumption that the spectral 
energy distribution results from continuous star formation 
and dust extinction. This assumption is justified by the fact 
that the luminosity d ensity is averaged ove r a large volume. 
According to models (iLeitherer et alll999l) . the unreddened 
s.e.d. can be reasonably approximated in the range 1200 A - 
2000 A by a power-law of slope a = —0.1 with current values 
of IMF and metallicity. Such a slope (—2.1 for s.e.d. per unit 
wavelength ) is observed in sta r-forming galaxies with low 
extinction iCalzetti et al.lll994l) . Models also show that this 
slope can be extrapolated shortward of 120 A but recent 
observations in the range 900 A - 1200 A iLeitherer et alJ 



4 LYMAN DISCONTINUITY 

A possible constraint on the parameter a is the need for the 
luminosity density at 900 A derived from the observed Ha 
luminosity density to be compatible with that observed at 
2000 A . The Lyman discontinuity predicted by stellar pop- 
ulation models is a key element in such a comparison. As it 
is cur rent practice llBruzual fc Charlotll993l : ILeitherer et alJ 
1999), the Lyman discontinuity is evaluated between 1000 
A and 900 A in order to avoid the influence of the severe 
line-blanketing by hydrogen towards the series limit at 912 
A. The luminosity density a t 1000 A is di r ectly related to 
that observed at 2000A bv lSullivan et all (l2000l) . through 
the intrinsic slope of —2.1 (per wavelength unit) and the 
extinction adopted A2000- Under current conditions, the Ha 
luminosity density o bserved and corrected for extinction by 
iGalleeo et aD <1995h (taken as 1.76 x 10 39 ergs s -1 Mpc -3 
for H = 70 km s -1 Mpc -1 ) implies a density of 1.3 x 10 51 
s -1 Mpc -3 ionizing photons, and in turn a luminosity den- 
sity of 0.67 x 10 38 ergs s -1 A -1 Mpc -3 at 900 A (assuming 
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a s.e.d. flat per wavelength unit in the Lyman continuum 
with no flux below the Hell edge). 

The resulting values of the Lyman discontinuity are 
listed in Table 1 as a function of A2000 • Most of these values 
are larger than the amplitude of 4 predicted for continu- 
ous star formation by curren t evolutionary synthes is mod- 
els jBruzual fc Chariot! Il993l : iLeitherer et alJll999ft . They 
would match this prediction for unrealistically low value of 
the extinction and in a domain of Fig. 1 where the integrated 
galaxy light would be larger than observed. 

However, a number of possibilities have been identified 
that, alone or in combination, would better match the eval- 
uation of the Lyman discontinuity with the amplitude pre- 
dicted by models. 

(i) If a fraction / of the ionizing photons escape or , more 
likely , are trapped by dust before ionization (e.g. Lnoui 

lUln ouc. Hiras hita fc KamavdbOOlllChariot et alJboO, 

DOT3ite^tHuT2003fl . the evaluation of the Lyman discontinu- 
ity would be multiplied by /. It would match the predicted 
value within the permitted range of background radiation 
from galaxies with a plausible value / = 0.7. 

(ii) The evaluation of the Lyman d iscontinuity would be 
decreased if the value adopted from iGallego et alJ <ll995h 
underestimates the Ha lu minosity density. Arguments have 



Table 1. Values of the Lyman discontinuity resulting from ob- 
servations as a function of the adopted A2000 extinction. 



been given in that sense bv lJones fc Bland-Hawthorn 



and higher values ha ve been reported Jc^^^^ _^ 
iGlazebrook et ai"ll2003t) : recently the value of TG^fleTOet^ll 
jl99a) has been revised upward by a factor of 1.6 
iPerez- Gonzalez et alJl2003f) . The evaluation of the Lyman 
discontinuity would be divided by the same factor and would 
be brought closer to the prediction. 

(iii) The possibility of an increase of the value predicted 
by models has also be examined. The amplitude of the stel- 
lar Lyman discontinuity is rather independe nt of metallicit y 
and only decreases at very low metallicity llSchaererl["2 003). 
Recent models including non- LTE effects <Schaererl Tl998: 
ISmith. Norris fc Crowtheil l2003fi tend to produce smaller 
discontinuities. In contrast, a relative content in massive 
stars reduced wit h respect to the curren t IMF would increase 
the discontinuity iLeitherer et al Il999l) . By far, a departure 
from the assumption of a continuous star formation would 
be the most radical way for increasing the stellar Lyman dis- 
continuity; we have no possibility to exclude an aging burst 
contribution in the composite population that would make 
the stellar Lyman discontinuity larger and consistent with 
the evaluation. 



5 DISCUSSION AND DUST EXTINCTION 

If the constraint from the Lyman discontinuity is ignored 
or assumed to be solved by the effects listed above, the 
predictions of the integrated light from galaxies are com- 
patible with the observations for a large range, but rela- 
tively low values, of the A2000 extinction (Fig. 1). Most of 
these values are lowe r than the extinction of 1.3 reported by 
ISullivan et al.l i2000t) as affecting their luminosity density. 
Taken at face value, the extinction A2000 = 1-3 would imply 
an evolution rate 7 ^ 3.5; a significant adjustment for the 
Lyman discontinuity and the lower bound of the UV inte- 
grated light would be required if 7 ~ 3.5 — 4. Such a definite 
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Values of the Lyman discontinuity are dependent on the evolution 
rate because the UV luminosity density had to be scaled back to 
z ~ (with the factor 1/(1.15)"') for comparison with the Ha 
luminosity density. 
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Figure 1. Predictions of the integrated light of galaxies at 1595 
A in units of photons cm~ 2 s" 1 A -1 ster" 1 as a function of the 
extinction A2000 f° r three values of the evolution rate, 7 = 1.5, 
2.5 and 4 from bottom to top (Ho = 70 km s —1 Mpc -1 , S1m=0.3 , 
=0.7). The range of values measured by G ardner et alj (2000'l 
is indicated as the two horizontal dashed lines. 



conclusion (7 ^ 3.5), however, is mitigated by a number of 
arguments. 

(i) The A2000 = 1-3 extinction compares well with 
the average UV extinction estimated in starburst galax- 
ies by a number of authors (e . g. iBuat fc BurgarelL 
I l998t iMeurer. Heckman fc Calzettil ll99Sl ICalzettil 1200 
iBell fc Kennicuttl I2001IL When it comes to UV selected 
samples of nearby galaxies, lower average extinction has 
been reported, especially when this extinction is eval- 
uated from the dust emission of the galaxies. UV se- 
lected samples are expected to contain normal, quies- 
cent star-forming galaxies in addition to starburst galax- 
ies as should be the case for the UV luminosity den- 
sity discussed here. Th e average UV ex t inctions reported 
range from 0.6 to - 1 <Buat fc Xulll996l: IBuat et al.lll999l : 
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llglesias-Paramo. Buat fe Donasl2o'o3h . Such values would be 
compatible with evolution rate 7 < 3. 5. 

(ii) The extinction of 1.3 found bv lSullivan et alJ |2000) 
for their luminosity density is surprisingly close to the aver- 
age performed on the individual extinctions, whereas a lower 
value is expected from a luminosity- weighted average. 

(iii) In addition to the A2000 extinction itself and the issue 
of th e universality of the starburst obscuration curve feelll 
120021) . the ratio fc(1000) /fc(2000) is a factor in the calculation 
of the integrated light from galaxies (cf equation 5). Given 
the uncertainties on the absorption curve, especially below 
1200 A, a value lower than the adopted 1.565 cannot be 
ruled out; this would raise the series of curves in Fig. 1 in 
the domain of large extinction and would make again the 
data compatible with a lower evolution rate. 



6 CONCLUSION 

We have tried to put together three quantities deduced from 
independent measurements, the accumulated far-UV light of 
galaxies (at ~ 1595 A), the luminosity density of galaxies in 
the far-UV (at 2000 A) and at Ha. As the star formation is 
likely to be continuous over large volume, we assume that the 
cosmic spectrum, i.e. the luminosity-scaled spectra summed 
over all galaxies, can be written in the far-UV as a simple 
power-law with an exponent depending only on the dust ex- 
tinction. In these conditions, the integrated light from galax- 
ies is predicted as a function of the cosmic evolution of the 
star formation rate and average dust extinction. This ap- 
proach will offer the possibility to pin down the cosmic evo- 
lution rate when the UV integrated light from galaxies is 
better known and the dust extinction better constrained in 
the UV (as with the GALEX survey). 

Data available at the moment imply an evolution rate 
at low z steeper than (1 + z) 3,5 . A shallower evolution is 
possible if the properties of the luminosity-weighted dust 
extinction are less extreme in the UV than those based on 
starburst galaxies. The latter possibility, with an average 
extinction ^2000 < 1-2 and/or a red dening law shallower a t 
short wavelengths than predicted bv lLeitherer et alJ 12002ft . 
cannot be ruled out and is supported by recent trends found 
in normal galaxies. 

The comparison between the luminosity density of 
galaxies in the far-UV (at 2000 A) and at Ha is compatible 
with the stellar Lyman discontinuity predicted by evolution- 
ary synthesis models if a fraction of the ionizing photons are 
trapped b y dust before ionizati on and/or the Ha luminosity 
density of lGallego et alJ Jl995T) is underestimated. 
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